
#include <iostream>
#include "fun.h"

int main()
{
    std::cout << "Test cuda Inputwise" << std::endl;
    {
    
    int size_x = 5;
    int size_y = 4;
    int bin = 5;
    hs::raster::RasterData<float> img(size_x, size_y);
    for (int i = 0; i < size_x * size_y; ++i) {
        img(i, 0) = (float)( i % bin );
    }

    printf("img: \n");
    printmat(img.view());

    hs::raster::RasterData<float> hist;
    histStatis(img, bin, hist);

    printf("hist:\n");
    printmat(hist.view());

    }

    std::cout << "Test cuda InputwiseIndex" << std::endl;
    {
    
    int size_x = 5;
    int size_y = 4;
    hs::raster::RasterData<float> img(size_x, size_y);
    for (int i = 0; i < size_x * size_y; ++i) {
        img(i, 0) = (float)(i + 1);
    }

    printf("img: \n");
    printmat(img.view());

    hs::raster::RasterData<float> dst_img;
    paddingZero(img, dst_img);

    printf("dst_img:\n");
    printmat(dst_img.view());

    }

    std::cout << "Done!" << std::endl;
    return 0;
}
